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We derive exact quantum expressions for the localization length L c for weak disorder in two- 
\ and three chain tight-binding systems coupled by random nearest-neighbour interchain hopping 

terms and including random energies of the atomic sites. These quasi- ID systems are the two- 
and three channel versions of Dorokhov's model of localization in a wire of N periodically arranged 
atomic chains. We find that L c 1 = 7V.£ 1 for the considered systems with TV = (1,2,3), where £ 
is Thouless' quantum expression for the inverse localization length in a single ID Anderson chain, 
for weak disorder. The inverse localization length is defined from the exponential decay of the 
two-probe Landauer conductance, which is determined from an earlier transfer matrix solution of 
the Schrodinger equation in a Bloch basis. Our exact expressions above differ qualitatively from 
Dorokhov's localization length identified as the length scaling parameter in his scaling description 
of the distribution of the participation ratio. For N = 3 we also discuss the case where the coupled 
chains are arranged on a strip rather than periodically on a tube. From the transfer matrix treatment 
tIh . we a l so obtain reflection coefficients matrices which allow us to find mean free paths and to discuss 

^ ■ their relation to localization lengths in the two- and three channel systems. 
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I. INTRODUCTION 



The localization length (L c ) of the quantum states is a fundamental parameter in mesoscopic physics. In particular, 
for quasi-lD disordered systems (wires) of finite length L, it not only sets the scale beyond which the electron states 
i ^ i ■ are effectively localized but it also determines the domain 

t-H ! 

> ' I < L < L c , (1) 

in which the conductance ql displays classical Ohmic behaviour, oc corresponding to a diffusive metallic 

' ■ regime. Here i denotes the elastic mean free path and 
CO 

o : 

L C ~N£ , (2) 



I 



where N oc y A is the number of scattering channels in a wire of cross-sectional area A. The metallic domain (1) 
does not exist for ID chains (N — 1) and for real wires with N >> 1 it requires the resistance to be less than 
some relatively large threshold value. These fundamental results have first been established by Thouless [1] and are 
*^ ' reviewed in [2,31. 

a: 

On the other hand, the notion of scattering channels itself is important since it has permitted the generalization 
of the well-known scaling equation for the evolution of the distribution of resistance (conductance) as a function of 
length in a ID chain [4], in terms of the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation [5,6] for a distribution 
of scattering parameters related to the conductance in quasi-lD systems. The DMPK equation together with the 
5—1 ' numerous results derived from it has been reviewed in [7-9] . 

a ' 

_ _ i 

Despite the important role played by the localization length in multi-channel disordered systems, microscopic an- 
alytic studies of it have remained scarce. This is the more surprising as L c is actually an intrinsically microscopic 
quantum parameter of fundamental importance, as recalled above. Note also that a first principles derivation of a 
relation between the localization length and the mean free path such as (2) would require separate calculation of both 
quantities in a disordered atomistic multi-channel system. 



Some years ago Dorokhov [10] has discussed a solvable model of multi-channellocalization consisting of N random 
tight-binding chains coupled by random nearest-neighbour interchain hopping terms and having random site-energies. 



Dorokhov's aim was to relax the assumption of isotropy of scattering parameters underlying the derivation of the 
DMPK equation [5-7] by replacing it by the weaker assumption of equivalent scattering channels [8] . After a fairly 
sophisticated analysis, which we find difficult to follow, Dorokhov arrives at an evolution equation for the distribution 
of scattering variables (participation ratio [11]) which involves a single microscopically defined scaling parameter, 
which he identifies with the localization length [10]. Dorokhov's expression for the localization length for weak disor- 
der is independent of the number of channels [Equation (12) below], which seems surprising! On the other hand, the 
popular models of transport and localization in multichannel wires such as the Thouless tunnel-j unction model [1,2], 
the random matrix and maximum entropy models [5-8] and the non linear sigma model [9], do not address detailed 
(discrete) microscopic models with specified disorder. 

In a recent paper [12] hereafter referred to as I, the author has derived exact analytical expressions for localization 
lengths for two- and three chain tight-binding systems with random site energies but constant nearest-neighbour in- 
terchain (transverse) and intrachain (longitudinal) hopping parameters, for weak disorder. In this model the channels 
are generally non equivalent, being associated with distinct channel-wavenumbers in the absence of disorder [12]. The 
localization length is defined, as usual, by the rate of exponential decay of the conductance. The conductance is 
determined using a transfer matrix approach for obtaining the amplitude transmission coefficients entering into the 
multi-channel Landauer formula. 

Motivated by our doubts about the correctness of Dorokhov's result (which, in particular, is incorrect for a ID 
chain) , we have reconsidered the calculation of the localization length for the case of two- and three equivalent chan- 
nels in his model, using the exact transfer matrix method for weak disorder developed in I. In view of the importance of 
Dorokhov's miscoscopic model in the context of scaling theories for probability distributions of transport parameters 
in quasi- ID systems it seems important to dispose of an accurate independent description of the localization length. 
On the other hand, the related analysis of reflection matrices will allow us to calculate mean free paths for the two- 
and three channel systems in the Born approximation and thus to test Eq. (2). 

In Sect. II we recall the Schrodinger tight-binding equations for Dorokhov's model for the case of two- and three chain 
systems. Dorokhov's model corresponds to periodic boundary conditions for the chains i.e. it describes equidistant 
chains arranged parallel to the axis on a tube. We also consider an alternative three-chain model with the parallel 
chains arranged on a planar strip, which corresponds to using free boundary conditions for the chains which are 
now non equivalent. In Sect. Ill we summarize the main points of the determination of the transfer- and scattering 
matrices in these models. Using the same parameterization of the matrix elements as in I, we limit ourselves to the 
definitions of these parameters in terms of the random site energies and hopping rates in Dorokhov's model, referring 
to I for presentation of the explicit forms of the corresponding matrices. The results and some concluding remarks 
are discussed in Sect. IV. In particular, we allude to a recently studied [13] weakly disordered multichain model 
including both interchain and intrachain nearest-neighbour random hopping but no site energy disorder. This model 
generalizes a well-known ID random hopping tight-binding model in which a derealization transition has been found 
at the band centre [14]. We give an exact expression for the localization length in this ID model, which readily reveals 
the derealization transition in the middle of the energy band. 

II. MICROSCOPIC MULTI-CHANNEL MODEL 

The TV-chain Dorokhov model [10] of a wire consists of parallel linear chains of Nl disordered sites each (of spacing 
a = 1 and length L = N^a) connected at both ends to semi-infinite ideal (non-disordered) chains constituting the 
leads. The sites on a given chain with its associated non-disordered parts are labelled by integers 1 < m < Nl in 
the disordered region and by m < 1 and m > Nl in the non-disordered ones, respectively. The disordered chains are 
coupled to each other by random hopping rates (transverse hopping) with vanishing mean values and, correspondingly 
the non-disordered chains are decoupled. The system is described by the tight-binding Schrodinger equation 

N 

^i +1 +i>L 1 + J2 £ n^ J n = Er n ,i=l,2,... N,l<n< N L . (3.a) 

3=1 

=£V4,™< lorn> N L . (3.b) 

Here E is the energy and tp^ denotes the amplitude of the wavefunction at a site m on the ith chain; e" = e l n is the 
random energy at a site n on chain i while 



2 



p ij — p ji 



(4) 



is a random symmetric hopping parameter between a site n on chain i and the corresponding nearest-neighbour site 
n on chain j. The above energies, including E, are measured in units of a fixed nearest-neighbour matrix element for 
hopping along the individual chains (longitudinal hopping). The random site energies and hopping parameters are 
assumed to be identically distributed independent gaussian variables with vanishing mean and correlation [s^ = (eo) 2 ] 

( £ n e m) = £ $i,j$m,n (5. a) 

(4 = 4 S m , n {5ij,6j, q + 5i, q 5 jtP ) . (5.b) 

We note that Eqs. (3. a) describe a collection of coupled chains of fixed separation, a, arranged parallel to the axis on 
a tube, which corresponds to periodic boundary conditions (pbc) for the chains. In the absence of disorder the chains 
are independent and equivalent and (3. a) shows that they all couple in the same way to the disorder. Therefore these 
independent chains define N equivalent scattering channels [10]. 

We now specialize to the cases of two- and three chain systems which are the object of this paper. For N — 2 and 
N = 3 Eq. (3. a) may be written 



f^n+l + IPI-A _(E- 4 - £l \ 2 \ _ 

U + i + ti-i) { -4 1 E - el) U 2 J ' 



For completeness's sake, we also consider, for N = 3, the case where the parallel chains are arranged on a planar strip 
which corresponds to free boundary conditions (fbc). In this case the Schrodinger equation is 



I = I -4 1 E-ei -4 3 I I CI • (8) 




Clearly, in this case, the channels are non-equivalent, but nevertheless well-defined. 

As in I, we shall determine the inverse localization length from the rate of exponential decay of the conductance of 
the disordered wires [1,2,15], 



where averaging over the disorder may be used, as usual, because of the self-averaging property of In g. The conduc- 
tance is given by the Landauer two-probe conductance formula [2,3], 



g = ^Tr(it + ) , (10) 



where t is the transmission matrix 



in 

*21 



ti2 ... tlAfX 
tN2 ■ ■ ■ tNN J 



(11) 
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where tij denotes the amplitude transmitted in channel i at one end of the wire when there is an incident amplitude 
in channel j at the other end. 

We close this Section by recalling the result for the localization length obtained by Dorokhov [10] for an TV-channel 
wire described by (3. a). In the notation of (3. a) and (5.a,b) it reads 



L r = 



A — E 2 



(12) 



which is independent of N . This surprising result follows by combining the expression for the localization length 
obtained from the scaling equation for the distribution of the participation ratio in the first equality of (6.26) in 
Ref. [10], with the definitions (2.9), (2.8) and (2.2). The Eq. (12) will be discussed further in Sect. IV. 

III. SUMMARY OF ANALYSIS 

Here we summarize the analytic study of transfer and scattering matrices for weakly disordered two-and three 
channel systems of I, as applied to the model of Sect. II. From this we obtain explicit results for the various intra- 
and interchannel transmission and reflection coefficients which we shall use for finding the localization length (9) and, 
as a check of our results, for verifying explicitly the current conservation property. The choice of similar notations to 
those used in I will allow us, conveniently, to refer to I for the explicit forms of the above matrices. 

A. Transfer matrices 

Transfer matrices, Y n , for thin slices enclosing only a single site n per channel of the system described by (6-8) arc 
defined by rewriting these equations in the form 



\ \ J 



IP 2 



V : 7 



(13) 



where 



Yn = X[ 



On 



1 

-el 1 
V 



-1 



-12 





E-e 2 n 
1 



£ 12 



(14) 



for N = 2, and 



Y n = 



for N = 3, with e\l = e%, and 



( E 


-4 


-l 




~Xn 3 






i 





















E-el 


-1 -£ 23 













1 










31 
An 





p 32 


E-el 


-1 


V 











1 


0/ 




Y n = 


X 'n 


Y i3 = 
i An 


for fbc , 





(15) 



(16) 
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Y n = X'l , x « = e « for pbc 



(17) 



In order to study the reflection and transmission of plane waves by a disordered wire we must use a basis corre- 
sponding to waves propagating independently from left to right and from right to left in the absence of disorder. Such 
a basis is provided by the Bloch waves supported by the system (3.a,b) for vanishing disorder. The Bloch waves are 
the solutions for real fc of the eigenvalue equation 



( <± 



± 

1,± 



,±*fc 



/ 1pn,± ' 

V>n-1,± 



(18) 



where Y" = X o, X' , Xq denotes the transfer matrices (14) and (15-17) in the absence of disorder i.e. s n = e\ = el\ = 
ejf = 0. The wavenumbers fc are given by 



2cosk = E , (19) 
for energies restricted to the band — 2 < E < 2. For definiteness we choose < k < ir so that the eigenfunctions 



< ± ~e ±irafc , (20) 
correspond to Bloch waves travelling from left to right and from right to left, respectively. 

In the absence of disorder, the transfer matrices (14) and (15) are diagonalized in the basis of Bloch wave amplitudes 
(20). In transforming (13) to the Bloch wave basis and performing, in particular, the corresponding similarity 
transformation (defined by the matrix W of the eigenvectors in (18) [12]) of the transfer matrix Y n for N — 2 and 
N = 3, respectively, we use the same parameterization for Y n = W~ 1 Y n W as in Eqs (22) (N — 2) and (23) (N = 3) 
of I (where we now put k\ = k 2 — k% = k). For the models (6-8) the parameters introduced in I are found to be given 
by: 



.12 



din = 



2sinfc 



«2r. 



2sinfc 



2sinfc 



(21) 



for N = 2, 



air 



2sinfc 



2 sin k 



hn = 



2 sin k 



fn 



,12 



2 sinfc 



-23 



2 sinfc 



(22) 



for N = 3 with fbc, 



n n t n 

" 2 sin k ' 3 ™ 2 sin fc ' 2n 2 sin k 



23 



fn — r, . , ? 9n — Pn — . , t d n — q n — . , (23) 

2 sin fc 2 sin fc 2 sin fc 



for N = 3 with pbc. 

Finally, the transfer matrices of the disordered wires of length L = N^a are the products of Bloch wave transfer 
matrices associated with the Nl individual thin slices n, 



5 



Nr. 



Y, 



= n 



n=l 



(24) 



For weak disorder it is sufficient to explicitate (24) to linear order in the random energies e l n and e 1 ^ for the purpose 
of studying averages to lowest order in the correlations (5.a,b). These correlations imply indeed that different slices 
in (24) are uncorrelated. The Bloch wave transfer matrices are given explicitly by Eqs (30) (N = 2) and (32) (N = 3) 
of I, with the parameters defined in (21-23) above and the wavenumbers ki,k 2 , fc 3 replaced by k in (19). 



B. Scattering matrices 



The scattering of plane waves (reflection and transmission) at and between the two ends of the random quasi- ID 
systems is governed by the S'-matrix, 



f++ 



where 



(25) 



*11 ^12 



(26) 



and 



P±=F 



'11 '12 



'21 



' 22 



(27) 



Here tf + (t~~) and r ij + ( r tj~) denote the transmitted and reflected amplitudes in channel i when there is a unit flux 
incident from the left (right) in channel j. Left to right- and right to left directions are labelled + and -, respectively. 
The S'-matrix expresses outgoing wave amplitudes in terms of ingoing ones on either side of the quasi-lD disordered 
wire via the scattering relations 




Here / and /' (0 and 0') denote ingoing (outgoing) amplitudes at the left and right sides of the disordered region, 
respectively. It follows from current conservation that e.g. for a unit flux which is incident from the right in channel 
i one has 



Ed *7rl 2 + I rj+ | 2 ) = 1 . (29) 

Likewise, one has also 



N 



Di*: 

i=i 



++ 

ji 



+ r 



= 1 



(29.a) 
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As shown in I, the wavefunction amplitudes at sites n and n — 1 in the Bloch representation correspond to waves 
travelling from left to right and from right to left, respectively. Like in I, we thus rename the wavefunction amplitudes 
in the Bloch representation of the transfer equations (13), 



w- 1 



/<4-l\ 

v ; / 



YnW- 1 



V : J 



by defining e.g. 



W~ 



f a tn-l\ 
a l.n-l 



(30) 



(31) 



Likewise, using a similar notation for wave amplitudes transferred from n = to n = Nl across the disordered wire 
of length L = N^a, we write the corresponding wave transfer equation in the Bloch representation, which follows by 
iterating (30), in the form 



(<l\ 






a t,L 


= Yl 


a h0 
a 2,0 






a 2,0 


V : ) 




V '■ J 



(32) 



The components of the out- and ingoing waves column vectors in (28) are thus a 1Q ,a 



2,0' • ' 



a i,o> a t,o> 



■ a N,0' a 



1,L> "2,1,' ' 



. . a 



N,L 



JV,0' a l,L' "2,L 



and 

■ a w l ' respectively. With the so defined vectors of outgoing and incoming amplitudes, 
the S-matrix is obtained by rearranging the equation (32) so as to bring them in the form (28). The details of this 
somewhat lengthy calculation are explicitatcd in I. The explicit forms of the scattering matrices, Eqs (46-47) and 
(48,48.a-48.f) of I, for N = 2 and N = 3, respectively, are expressed in terms of transfer matrix elements which are 
themselves defined in terms of general parameters given by (21-23) above in the case of Dorokhov's model. These 
S'-matrices readily yield the transmission and reflection submatrices in (25). 



IV. RESULTS AND CONCLUDING REMARKS 



The explicit expressions of the transmission- and reflection coefficients, \t^~\ and + | , in terms of the general 
parameters defining the transfer matrices in I are given in an appendix in I. By inserting the present parameter values 
(22) and (23) for the two- and three-channel Dorokhov models in these expressions and averaging over the disorder, 
using (5. a) and (5.b), we obtain the following results, exact to order e$: 



|irrl 2 ) = (fe-| 2 ) = i 



ZN L el 
4 sin 2 k 



(33) 



(|tr 2 -| 2 ) = (|t 2 Ti 2 > = ^ 



4 sin 2 k 



(34) 



l^ + | 2 > = 



4 sin 2 A: ' 



i,j = (1,2) , 



(35) 
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for TV = 2, 



(36) 



4 sin k 



(37) 



(\t- 12 -\ 2 ) = {\t- 21 -\ 2 ) = ^4- h , 

4 sm k 



(38) 



(fe-| 2 > = (K 3 Tl 2 ) 



4 sin^ fc 



2 7, ' 



(fe"l 2 ) = (I^Tl 2 > = o 



(39) 
(40) 



Ti + I 2 > - (K/l 2 ) = (k33 + l 2 > = (kr 2 + l 2 > - (l^+l 2 ) 



for TV — 3 with free b.c, 



(|rT3 + | 2 ) = (|r 31 + | 2 )=0 



(l^3 + | 2 ) = (l^| 2 ) = S| , (41) 



(42) 



(l^f| 2 ) = (&| 2 > = (|i3Tl 2 > = l 



5N L e 2 
4 sin 2 k 



(43) 



<l*«-| 9 > = 4^ ' W > M = (1,2,3) 



(44) 



|r,7T> 



A^e 2 , 
4 sin 2 fc 



, *,i = (1,2,3) , 



(45) 



for TV = 3 with periodic b.c. 

One readily verifies, as a check of the explicit results (33-45), that in all cases (TV = 2, TV = 3 with fbc and TV = 3 
with pbc) the current conservation property (29) is obeyed. 

Next, by evaluating the averaged traces (Tr{i~~(i~~) + ]), successively for the three models using (33-34), (36-40) 
and (43-44), respectively, we get: 



(Tr[t— (*—)+]) = 2 - , V = 2 



sin 2 k 

= - '""o° , N = 3 with fbc , 



7TV £ e 2 
4 sin 2 k 



9TV L e 2 



4sin 2 fc ' 

For the inverse localization lengths defined in (9-10) we then obtain: 



TV = 3 with pbc . 



(46) 
(47) 
(48) 
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F 2 



L c 4 sin k 
1 _ lei 



L c 24 sin 2 k 
1 3e 2 



L c 8 sin k 



, for N = 2 , (49) 
, for N = 3 with fbc , (50) 
, for N = 3 with pbc . (51) 



These expressions are exact to order el for weak disorder. 

It is instructive to compare (49-51) with the localization length, £, for weak disorder in a one-dimensional chain with 
random site energies. In this case Thouless [16] obtained the exact expression 

7 = ^1 ' (52) 
4 8 sin k 

which has been rederived in I (see also [17]) using transfer matrices. We observe that the inverse localization lengths 
for pbc in (49) and (51) take the values 2/£ and 3/£ for N — 2 and for N = 3, respectively. The constant value 
(12) obtained by Dorokhov [10] for arbitrary N differs qualitatively from these exact results. Using (19) and (52) 
Dorokhov's expression may be written 



1 4 

— = — , N arbitrary . (53) 

In fact, the exact expressions (49) and (51), together with the ID expression (52), suggest that the actual form for 
the inverse localization length for weak disorder for arbitrary N could be 

We also note, incidently, that in our analysis of the two- and three-channel Dorokhov models the localization length 
reduces precisely to the ID result (52) in the limit of no interchain hopping (e% = 0), as expected. Indeed, for e l J = 0, 
we have t^-~ = r^- + = for i ^ j and from the explicit expressions of the random transmission coefficients in the 
appendix of I we get, using (21-23) and (5. a), 



(Tr[i-(i-)+])=2-^4r,N = 2 
2 sin k 



= 3(l-^4-),N = 3 , 
V 4sin 2 /cy 

which both lead to l/L c — l/£. 

On the other hand, the above results for reflection coefficients may be used for obtaining explicit expressions for 
mean free paths in the few-channel quasi-lD systems. The mean free path for an iV-channel wire is defined by [8,18] 



£-ssf£<M'> - (55) 

We then obtain successively from Eqs. (35), (41-42) and (45) 

(56) 

(fbc) , (57) 
(pbc) . (58) 
In the one-dimensional case one gets similarly, by returning from (52) to the determination of the reflection coefficient, 



f 


el 




2 sin 2 k 


1 


7 e 2 
i e 




12 sin 2 k 


1 


3 e 2 


T 3 


4 sin 2 k 
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The expressions (56-59) correspond to the Born approximation of impurity scattering. By comparing (56-59) succes- 
sively with the localization lengths in (49-52) we find that in all cases 



L C = 2£ N , TV = 1,2, 3 . (60) 

We also note that a similar calculation of mean free paths for the two- and three channel wire models with constant 
interchain hopping rates discussed in I [12] also leads to Eq. (60) for N = 2 and N = 3. The localization lengths for the 
multichannel systems in I are given by Eqs (58), (73) and (86) of that reference, respectively, and the corresponding 
reflection coefficients entering in (55) above, are given by Eqs (52-57), (67-72) and (83-85) of [12]. 

The Eq.(60) for the one-dimensional case coincides with the relation between the localization length and the 
mean free path derived by Thouless from kinetic theory [19]. Our treatment thus establishes a similar exact re- 
lationship for two- and three-channel systems both for Dorokhov's model and for the model with constant in- 
terchain hopping in I. The exact universal expression (60) differs qualitatively from Eq. (2) discussed earlier, 
mainly for N >> 1 [1,2,6,7,9,10], and does not suggest the existence of a well-defined diffusive (metallic) regime, 
In « L « L c , in few-channel systems. We recall that in the above references the mean free path is introduced 
as a fixed length scale beyond which metallic diffusion takes place (when it is not inhibited by localization). Our 
microscopic analysis yields explicit expressions both for localization lengths and for mean free paths. 

Finally, from (50) and (51) it follows that the difference in transverse boundary conditions for the corresponding 
three-channel models has only a minor influence on the localization lengths. 

The transfer matrix approach discussed in I may also be applied for studying the derealization transition which has 
recently been found at the band centre in weakly disordered multi-chain systems including both nearest-neighbour 
inter- and intrachain random hopping terms but no site energy-disorder [13]. This derealization transition exists 
already in a one-dimensional chain with random hopping, as has been known for some time [14]. In this case it may 
be readily revealed by studying the localization length L^ 1 = — limL^ 00 (2L) _1 (ln |i ±:t | 2 ) of the chain. Consider the 
Schrodinger equation 



(1 + r) n )(tp n+1 + tp n -i) = Eip n , (61) 

where rj n is a gaussian random nearest-neighbour hopping parameter (with (rj m rj n ) — "q^8 m>n ) measured in units of 
the non random hopping parameter. From a transfer matrix analysis of (61) similar to that used for obtaining the 
transmission coefficient and the corresponding localization length (52) for a weakly disordered Anderson chain [12,17] 
we get 



1 = vjc^k E = 2cQsk (g2) 
Li c I sin k 

This expression, which is exact to order ry 2 , displays the divergence of the localization length in the middle of the 
energy band, E = 0. 
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